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Abstract 

We examine the role of resonant coupled channel final state interactions (FSI), as well as weak 
annihilation and exchange contributions in explaining all the two body hadronic D —>• PP decay 
modes data. In the un-unitarized amplitudes we include modified Wilson coefficients with non- 
factorizable corrections as parameters. For the hadronic form factors, the z-series expansion method 
is used to get the q 2 dependence. The FSI effects are incorporated via a phenomenological approach 
with widths of resonances to various channels taken from observations where available, and others 
as additional parameters to be determined from fits of all the theoretical rates to the measured ones. 
Our results for the rather hard to explain D° —>• K + K~ are in agreement with measured 

values. We demonstrate that both weak exchange as well as FSI effects are required to get the 
correct branching ratio for the D° —> K°K° mode. Using our unitarized amplitudes we evaluate 
the strong phase difference between the amplitudes for D° —* K~tt + and D° —> K + tt~ and find it 
to be in complete agreement with the recent BES III result. 
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I. INTRODUCTION 


With the advent of the B factories, along with even the Tevatron having sufficient lumi¬ 
nosities to perform excellent heavy flavour studies and more recently with the data pouring 
in from LHCb and BESIII, quark flavour physics is changing its role from being the ground 
for confirming the Kobayashi-Maskawa mechanism for CP violation, to providing the ob¬ 
servables that can test or constrain physics beyond the Standard Model. With improved 
precision measurements the goals of hadronic weak decays will cease to be, from seeking 
observables free from hadronic uncertainties, to actually precisely understanding the QCD 
effects. 

Along with the plethora of data on semi-leptonic, hadronic, rare decays, CP asymmetries 
and mixing in B decays, discovery of charm mixing and hints of CP violation in the charm 
sector, resulted in non-lcptonic charm meson decays being a focus of attention, in the last 
few years [1, 2], Since CP violation in charm, direct as well as in mixing, is expected to 
be negligible, any hint of CP violation in charmed mesons is expected to be from physics 
beyond the standard model and hence charm may be instrumental in providing insights into 
new physics [3-11]. In fact, earlier FOCUS [12], CLEO [13] and various other collaborations 
had produced many interesting results in the charm sector. The 3.2 a hint [14, 15] of a 
difference of CP asymmetries between the singly suppressed K + K~ and charmed 

decay modes resulted in a large volume of work [16-32], mostly using different models of 
New Physics, to explain the result. The hint has since then been slowly moving towards 
zero and currently there seems to be no evidence for any direct CP violation in charm in 
any mode [33]. Even though this hint for the CP asymmetry slowly disappeared, from all 
the recent work done on the charm decays to two pesudoscalar modes motivated by this 
hint, it was clear that it is critical to first understand the observed branching ratios of all 
the charmed hadronic decay modes well, within the Standard Model, before any observation 
of an anomalous rate or any new CP asymmetry can be claimed to be, due to the presence 
of New Physics. 

However, this is not an easy task. The mass of the charm quark (1.275 GeV) makes 
it very difficult to come up with a proper theoretical technique for calculation of hadronic 
charmed meson decays. The charm, unlike the bottom quark, is not sufficiently heavy to 
allow realization of the infinitely heavy quark limit. Therefore, the well-known theoretical 
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approaches based on QCD, for example, heavy quark effective theory [34, 35], QCD fac¬ 
torization [36, 37], the perturbative QCD approach [38-41] and the soft-collinear effective 
theory [42], which lead to very satisfactory predictions for B decays, cannot be used to 
explain data in the case of charmed mesons. Furthermore, the charm quark is also not light 
enough for a chiral expansion to be applicable. 

In the absence of any other reliable and effective theoretical methods, the factorization 
approach is still one of the most successful ways to study two-body charm meson decays [43, 
44], However, it is well known now that in the naive factorization approach, calculation 
of Wilson coefficients of effective operators faces the problem of y 5 - and renormalization 
scheme dependence. These difficulties can be overcome in the frame-work of the ‘generalized 
factorization approach’ where Wilson coefficients are effective and include important non- 
factorizable (NF) corrections [45, 46]. 

In the past there was another attempt to explain hadronic D decays using the so called 
large 1/N C (where N c is number of colour degrees of freedom) approach [47]. It was observed 
that dropping Fierz transformed terms characterized by l/iV c can narrow the gap between 
predictions and observations up to a satisfactory level. The calculations based on QCD sum 
rules showed that Fierz terms were certainly compensated by the NF corrections [48-50]. 

There exists another model independent so called ‘quark diagram’ or ‘topological diagram’ 
approach in the literature [51-59] where all two-body non-leptonic weak decays of heavy 
mesons are expressed in terms of distinct quark diagrams, depending on the topologies of 
weak interactions, including all strong interaction effects. It is based on SU( 3) symmetry, 
and allows extraction of the quark diagram amplitudes by fitting against experimental data. 
However SU(3) breaking effects in charmed meson decays have been shown to be important 
and need to be carefully incorporated [60, 61] . 

The importance of final state interactions (FSI) in nonleptonic charm decays had been 
realized and discussed in several papers [62-65] in the early 80’s, where the authors had been 
intrigued by the anomalies in the observed branching ratios of the Cabibbo favoured (CF), 
neutral versus the charged Kit modes, the differing rates of the singly Cabibbo suppressed 
(SCS) K + K~ and 7r + 7r - modes, followed by measurements of rates of few other modes that 
had unexpected suppression/enhancement. Many of these were conjectured to be due to 
FSI. Surprisingly, even in the last couple of years, in many of the papers that worried about 
the charm CP asymmetry problem, these old puzzles were still considered unresolved. 
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Even for the case of hadronic B meson decays, the role of FSI’s is being examined rather 
carefully in the last few years [66, 67]. The mass of the charmed meson lies right in the 
heart of the resonance region. Hence, resonant final state rescattering is bound to play a 
bigger role in the two body hadronic charm decays and needs to be evaluated. Of course 
dynamical calculations of these long distance effects are not possible and hence they can 
only be determined phenomenologically after comparison of the theoretical estimates with 
experimental data. Unitarity constraints play an important role in providing the theoretical 
estimates. 

Another contribution in hadronic two body decays that has been debated for a long 
period is that of the weak annihilation and exchange diagrams. One of Rosen’s proposal [68] 
had been that the W exchange diagrams may be large and since this appears only in D° 
and not in D + decays, it could account for the difference in the lifetimes of these two 
mesons. Bigi and Fukugita [69] had then proposed several D and B meson decay modes 
that could be the smoking gun signals of the W-exchange contributions and yet, when the 
mode D° —>■ (j)K° was observed, it was argued [70] that it could have been generated from the 
decay mode D° —> K*r /, with this final state rescattering to the cpK° mode. Annihilation 
type contributions along with FSI’s were incorporated in the hadronic two body vector- 
pseudoscalar modes of charmed meson decays in Ref. [71]. Studies using the quark diagram 
approach of Ref. [72] had also indicated that annihilation type contributions are needed to 
explain the observed data. 

In this paper, we study the role of FSI in the two body D (.D here can be any of the 
D°, D + , or Df) meson decays. We assume that FSI effects are dominated by resonance 
states close to the mass of D mesons. In fact, there exist isospin 0, 1 and 1/2 resonances 
near the D mass, that may contribute to rescatterings among different channels in these 
respective isospin states and enhance/suppress some of the decay rates. In the next section, 
we give the formalism for the calculation of the un-unitarized amplitudes, using a modified 
factorization approach, where the effective Wilson coefficients include NF corrections. This 
is in analogy with the QCD factorization approach of Beneke-Neubert for hadronic B meson 
decays [73], where however, using the hard scattering approach, the NF corrections are 
calculable in heavy quark approximation. However, for charm, since this approximation 
fails, these NF corrections are not calculable and are left as parameters. We also indicate 
our parametrization of the annihilation contributions and discuss our inputs: the decay 


4 



constants, and the form factors, for which we have used a z-series expansion approach. In 
Sec. Ill, we discuss the need to incorporate additional long distance FSI effects and show 
how this can be done with a K matrix formalism for coupled channels. Using the observed 
widths, masses and known decay rates of the resonances to the various channels to evaluate 
the diagonal elements and leaving the unknown elements of the K matrix as parameters, 
the unitarized amplitudes are calculated (as discussed in [71]) to estimate the branching 
ratios of all the SCS, CF and doubly Cabibbo suppressed (DCS) D —> PP decay modes. In 
Sec. IV, we list the isospin decomposition of all the decay modes, the parameters that need 
to be determined from our fits as well as the errors in theoretical inputs used. We list all the 
branching ratios after our numerical \ 2 fits as well as the values of the fitted parameters. 
Finally we conclude in Sec. V. 

II. THE UN-UNITARIZED AMPLITUDES 

A. Weak Hamiltonian and Wilson coefficients 

The study of weak decays of charmed mesons to two body hadronic modes necessarily 
requires a careful evaluation of the strong interaction corrections. The weak effective Hamil¬ 
tonian may be expressed in terms of coefficient functions, which incorporate the strong 
interaction effects above the scale /x ~ m c and the current-current operators as: 

Uw = + C 2 {n)0 2 {n)} + h.c. . ( 1 ) 

where, Gp is the Fermi coupling constant, C\ and C 2 are the Wilson coefficients and the 
operators are, 

0\ = ( U a q 2 a)v-A((hpCp)v-A 
O 2 = ( U a q2/3)v-A(<h/3C a )v-A ■ 

a and (3 in the above are colour indices, while q \, q 2 can be either the d or the s quark. 
The quark diagrams dominantly contributing to the branching ratios of D —> P 1 P 2 [59] are 
the colour-favoured Tree amplitude T, the Colour-suppressed amplitude C, the W-exchange 
amplitude E and the W annihilation amplitude A, shown in Fig. 1. 

Penguin contributions in charmed meson decays are highly suppressed as the dominant 
down type quark contribution to the flavour changing neutral current c —> u transition is 
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from the b quark which is accompanied by the presence of the tiny product, V* b V u b of the 
CKM matrix elements. Hence, the two operators in Eq. (1) are sufficient for calculating the 
amplitudes and branching ratios of the D —> PP modes. 

In the naive factorization approach, the matrix element of the four-fermion operator in 
the heavy quark decay is replaced by a product of two currents. The amplitudes for the 
non-leptonic 2 body decay modes are then the product of a transition form factor and a 
decay constant. However, NF corrections must exist; while such corrections for scales larger 
than yU are taken into consideration in the effective weak Hamiltonian, those below this scale 
also need to be carefully incorporated. In the QCD factorization approach for B meson 
decays [36, 37, 74], these NF corrections are handled using the hard scattering approach, 
where the vertex corrections and the hard spectator interactions are added at the next to 
leading order in a s and its accuracy is limited only by the corrections to the heavy quark 
limit. But, in the case of charm decays, where the heavy quark expansion is not a very good 
approximation, it is best to parametrize these NF corrections and then determine them by 
fitting the theoretical branching ratios with the experimental data. In the diagrammatic 
approach of Ref. [59] also, either the Wilson coefficients themselves or the NF corrections 
appearing in the Wilson coefficients are determined from fits to data. 

Hence, we write the scale dependent Wilson coefficients, modified to include the NF 
corrections which are parametrized by Xi and Xi with their respective phases and (j ) 2 as, 

cii(h) = CM + CM + a(/i)yie*^ 1 ^ (2) 

a 2 (^) = C 2 (/i) + Ci(fi) + a(/i)y 2 e*^ 2 ^ • (3) 

The dominant Tree and Colour amplitudes for D —> P\P 2) where P\ is the final meson which 
carries the spectator quark, while P 2 represents the meson emitted from the weak vertex (as 
depicted in Fig. 1), are then written as: 

T ( c ) = V C KMaM(a 2 (v))fp 2 (m 2 D - m 2 Pl )F 0 DPl (m 2 p 2 ) , (4) 

where, fp 2 is the P 2 meson decay constant and P 0 DPl (mpJ denotes the transition form factor 
for D —> Pi evaluated at m 2 p 2 . We follow the prescription of Ref. [75] and choose the scale 
[i to be the energy release in individual decay processes rather than fixed at m c . This 
scale which is dependent on the final state masses, allows for SU(3) breaking, additional 
to that coming from different decay constants and form factors. This scale is taken to be, 
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q q 

(a) Emission or Tree (T) 



q q 

(b) Colour suppressed (C) 




(c) W-exchange (E) (d) W-annihilation (A) 


FIG. 1. The dominant quark diagram amplitudes 


p = — r|), w l iere 5 r 2 = m p 2 / m z) and A is another free parameter. A, xi , X 2 , 4>i 

and 02 are taken to be universal for all the decay modes and are fitted from experimental 
data. 


B. Weak Annihilation Contributions 

For a long time, W-exchange and W-annihilation contributions used to be neglected due 
to the so-called helicity suppression. However, observation of many decay modes of charmed 
and bottom mesons, which are possible only via the annihilation or exchange diagrams have 
indicated that these contributions could be substantial. These short distance weak annihi¬ 
lation effects were hence included in the diagrammatic approach. In principle these could 
result from rescattering, even in the absence of annihilation and exchange processes. In fact, 
weak annihilation topologies were assumed to be induced by nearby resonances through 
FSI’s in Ref. [72], The authors in [72] as well as Ref. [76] use SU(3) to relate the couplings 
of the final state mesons with the resonances. This leads to the result that the long dis- 
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tance W-exchange contribution can be induced by a tree amplitude, while W-annihilation 
can be induced by a color suppressed internal W-emission. The resonant FSI modify the 
W-exchange and W-annihilation amplitudes but the T and C amplitudes are unaffected. 
We emphasize that the assumption of SU(3), plays an important role in these results. In 
a most general coupled channel formalism, all contributions in the various channels will be 
affected by the resonant FSI, as will be shown in Sec III. We parametrize the W-exchange 
and W-Annihilation contributions in the amplitudes by y EI ' A ' 1 and estimate them from phe¬ 
nomenological fits to data. The exchange(annihilation) amplitudes are hence written as: 


E q ,s(A q , s ) ~ -^ v CKMCi(fi)(C 2 (fJ.))xf!s A) -j£f.bfpj p. 


a/2 


AW 


( 5 ) 


Since the initial charmed meson is annihilated and both the final mesons are produced from 
the weak vertex, after the production of a quark-antiquark pair from a gluon, these ampli¬ 
tudes are a product of the decay constants of the initial D meson (f D ) and that of Pi(/pi) 
and P 2 {fp 2 )- Apart from this, the strengths of the exchange(annihilation) amplitudes Xq^ 
are assumed to be the same for all modes and the subscripts distinguish the contributions 
of the pair production of the light quark-antiquark from that of the strange pair. Since the 
annihilation and exchange contributions are necessarily non-faetorizable, they depend only 
on C i j2 rather than the modified coefficients aq^. 

We wish to emphasize that for the case of decays to two pseudoscalar mesons, it has been 
shown in Ref. [73] that in the annihilation contributions the quark-antiquark pair produc¬ 
tion happens with gluon emission from the initial state quark. Hence, this contribution is 
independent of the FSI effects that obviously involve the final state quarks and are discussed 
in Sec. III. Hence inclusion of both weak annihilation/exchange contributions as well as FSI 
will not amount to double counting, as pointed out in Ref. [77]. 

The scale of the Wilson coefficients for the exchange and annihilation amplitudes must de¬ 
pend on both the mass ratios, rq^ = mp 12 /mn and is taken to be, [i = s/Arnpl T^rf)(T^r§) . 



C. Non-perturbative Inputs: Form Factors and Decay Constants 


We start by specifying our convention for the different mesons involved in our analysis: 


7T+ = — U(1, 7T- = dU, 7T° = , 

K° = ds, K° = —sd, K + = us, K~ = sii, 
D° = cu, D + = —cd, Df = cs. 


In the D —y P\ transitions, the matrix element of the vector current is written in terms 
of the form factors F + and F 0 as, 


2 2 
'rri ^ r rrf ^ 

( Pi(p')\qYc\D(p )) = F + (q 2 ){p^ + - - - 


, g") + F 0 (g 2 ) m ^ 2 m ^ g", (6) 

q* q 

where q = p — p'. The matrix element for the production of the second meson P 2 , is given 

by, 

(P 2 id) I1 0 ) = ifp 2 q„ • (7) 


Hence, in the amplitude of the non-leptonic two pseudoscalar decay modes of charmed 
mesons involving the product of the two matrix elements specified in Eqs.(6) and (7), only 
the transition form factor Fq appears. Transition form factors can in principle be exper¬ 
imentally measured from the semi-leptonic decays, however, in the massless lepton limit, 
only the F + (q 2 ) contributes to the semi-leptonic amplitude distributions. However, the 
semi-leptonic information is still useful, since at zero momentum transfer the form factors 
obey the kinematic constraint Fq( 0) = F + (0). The q 2 dependence of the Fq on the other 
hand is accessible only with massive leptons in the scmileptonic decays or in lattice sim¬ 
ulations. Simple and modified pole models have been widely used to parametrize the q 2 
dependence of the form factors, but these have poor convergence properties. Recently the 
z-expansion [78, 79] has been introduced as a model independent parametrization of the q 2 
dependence of form factors over the entire kinematic range and has been shown to have 
improved convergence properties. In this approach, based on analyticity and unitarity, the 
form factors are expressed as a series expansion in powers of z n , where z is a non-linear 
function of q 2 , with an overall mutiplicative function accounting for the sub-threshold poles 
and branch cuts, 

inf 

F (t) = -7 77 , T a k(to)z(t, to) k ■ (8) 
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The series coefficients and prefactors can only be determined from fits to lattice or ex¬ 
perimental data. In fact, CLEO collaboration has determined these coefficients for the 


D —> n,K,r) form factors from the semileptonic decays but, in the massless lepton limit. 
Hence, for F 0 (q 2 ), we use lattice results to determine the first two coefficients. 

In the Becirevic and Kaidalov (BK) ansatz [80], 

fota 2 )"-^-, 1 (9) 

where, m* D is the mass of the vector meson with flavour cd or cs, depending on the transition 

being c —» d or c —» s respectively. F 0 (0) and fd are parameters to be fitted to experimental 

data and in fact, in Ref. [82], the Fermilab and Lattice MILC collaborations, have fitted 

these parameters to CLEO-c data. The normalization /(0) and shape parameter determine 

the physical observables describing the form factors at large recoil and are given by, 

j/n\ — p /n\ __ p /n\ 1 _ (Mg ~ M|) dF 0 

/(0) _ F+{ 0) F 0 (0), ^ dq2 W =o • (10) 

Using these input parameters we can determine the first two coefficients of the series expan¬ 
sion in Eq. (8) for F 0 (q 2 ). 

Few details regarding the z series form factor expansion can be found in Appendix B. If, 
the series is rapidly converging, even two coefficients may be sufficient to determine the q 2 
dependent form factor F 0 (q 2 ). Equating the normalization and slope obtained using Eq. (8) 
to that obtained from the lattice parameters (/(0),/3), which in turn had been obtained by 
fits to experimental data, we can obtain the z-expansion series (up to linear order) for all 


the form factors. 


In Fig. 2 we show the plots for our results for D —> n and D —> K where the lattice input 
parameters from Ref. [82] are used. We would like to point out that our D —> tt and D —> K 
form factor values at q 2 = 0 are in very good agreement with that given in Refs. [83, 84] 
which are the most precise published calculations for D —> nlv and D —> Klu form factors, 
according to Lattice Review [85]. Further, the shape of F 0 (q 2 ) for D —>■ -it, K that we obtain 
after the z-series expansion are consistent with that of Ref. [86]. Moreover, the /(0) values 
in Table I for D —» n and D —> K are also in agreement with CLEO results of Ref. [87]. 

Regarding D —>• rj and D —> q' form factors, since in these transitions, only the non- 
strange component q q is involved, hence one expects, F 0 (q 2 ) for D — > rj q ,rj' q ~ D —>■ tt. For 


1 It has been pointed out in Ref. [81] that Fq, but not F+, can be modeled by a single pole. 
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FIG. 2. The q 2 dependence of the scalar form factors. The plot on the left displays Fo(q 2 ) for 
D —>• 7r transition, while that on the right is for D —>• K transition. 

D —>■ 77 , CLEO has determined F + (0)V c d = 0.086±0.006±0.001 using the semileptonic decay 
mode D + —> r]e + u e [ 88 ]. Hence we use this value to estimate F 0 (0). However, for 77', while 
the first observation of the decay mode D + —> r]'e + v e has been reported by CLEO in the 
same paper, but the form factor was not determined in this case, and hence we approximate, 
the Eo(0) for 7 / to be the same as that for 77 . Further, since it has been shown[89] that the 
form factors and particularly their shape is insensitive to the spectator quark, the shape for 
both 77 and 77 ', is assumed to be the same as that for the D —>■ 7 r case. 

The D s —> 77 , 77 ' have been estimated with some lattice studies using the (/(0),/3) values 
from a recent exploratory paper by Bali et. al. Ref. [90]. However, these have larger 
uncertainties, since even the lowest pion mass used is still far from the physical mass. Hence, 
for D s —>- 77 , 77 ' we take the form factors to be similar to that of D —> K and for D s —>■ K 
we take them to be similar to D —> n. Note that in all these cases the masses of the 
final mesons for each of the respective decay process are used in obtaining their z-expansion 
coefficients, the approximations are used only for the input parameters, /(0) and j3. Due to 
this uncertainty, we have added an additional 3% theoretical error to these form factors. 

Turning now to the decay constants, for n and K mesons, the f^.K are taken from the 
Particle Data Group(PDG) [91]. For the 77 and 77 ', following the method described in [92], it 
is assumed that the decay constants in the quark flavour basis, follow the pattern of particle 


I 1 1 1 1 I 


D—>K 


1 _ 1 _ 1 _1_ 1 _ 1 _ 1 _ 1 _ 1 _ 1 _ 1 
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TABLE I. Best fit values of BK parameters for the scalar form factors 


Decay 

/(o) 

P 

D — > 7T 

0.64 ±0.03 ±0.06 

1.41 ±0.06 ±0.07 

D-+ K 

0.73 ±0.03 ±0.07 

1.31 ±0.07 ±0.13 


TABLE II. z-expansion coefficients obtained after using the BK parameters in Table I 


Decay 

ao 

ai 

D —)• 7T 

0.19 ±0.02 

-0.41 ± 0.05 

D -> K 

0.08 ± 0.01 

-0.32 ±0.03 

D q 

0.06 ± 0.004 

-0.27 ±0.02 


state mixing. The 77 and rf are expressed as linear combinations of the orthogonal flavour 
states, 


r] q = —=. (uu ± dd ), and r] s = ss . 

a/2 

The physical states rj and rf are related to these flavour states by, 


( 11 ) 



where, the 77 — rj 1 mixing angle denoted by 0 , represents the sum of the ideal mixing angle 
and the 77 — 7 / mixing angle ( 6 ) in the octet-singlet basis, 0 = 6 ± tan -1 y/2. Hence the 
decay constants (form factors) f q and f s (F 0q and F 0 J corresponding to that for r] q and rj s 
(D —> r] q and D —> r} s ) respectively, are given by: 


ff = fq cos 0 , 

ff, = ~fs sin 0 , 


= F 0q cos 0, 

F s — 
r o v ~ 

—Fo s sin 0 . 

fp = fq sin 0, 

fp = fs cos 0 . 

% 

= F 0q sin 0, 

F s 

r o v , 

= F 0s cos 0. 


Various ratios of decay rates having rf in the final state with respect to that with 77 , for 
example, T(J/0 — > r]' p)/Y{J /'if — > 77 p), comparison of cross-sections of scattering processes 
for 7 r~p —> rfn with that of n~p — > rjn etc., had been used for a phenomenological fit for 
the decay constants as well as the angle 0 in Ref. [92] and had been widely used. Recently, 
Babar with more accurate data on two photon widths of light pseudoscalar mesons, did a 
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combined analysis [93] along with CLEO data to yield a mixing angle and decay constants 
with reduced uncertainties: = 37.66 ± 0.70, fy- = 1.078 ± 0.044 and = 1.246 ± 0.087, 
which are used in this work. 

With the above inputs, the un-unitarized amplitudes for all the two-body pseudoscalar- 
pseudoscalar (PP) modes: SCS, CF and DCS may be written and are listed on the next 
page. For the decay modes involving r/ and rf , to distinguish the case in which r/ q is the P 2 
meson of eqn (4) from that where r] s is the one, the notation used is: 

c h = 

c i = V C KMa 2 (v)fs(ml - 

c {' q = < ^=VcKMa 2 {fi)f q (rn 2 D ~ m 2 Pl )F° Pl (m 2 v ,), 

c !y s = ^VcKMa 2 {^)fs(m 2 D - m 2 Pl )F PPl (m^). 


Similarly, to distinguish the cases where rj q or r] s is the P\ meson, which incidentally 
appears in both Tree(T) and Colour Suppressed(C) amplitudes (unlike for the case discussed 
above), we use the notation: 

c H°( T v°) = C ^= v CKMa2(li)(ai(ti))fp 2 (m 2 D - m 2 )F°’ l ‘(m 2 p J, 

C%(T?°) = ^jfcKMO. 'M(aM)fn« - 

C%(T%) = P r OT „a 2 ( / .)(a 1 (rt)/p a (mJ, - m^F^tm^). 

III. FINAL STATE INTERACTIONS 

Final state interaction effects are incorporated using unitarity relations, where the contri¬ 
bution to any channel is a result of sum over all possible hadronic intermediate states. Hence 
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SCS Decays 


A(D° - 
A(D° - 

A(D° - 

A(D° - 

A(D° - 

A(D° - 

A(D° - 
A(D° - 
A(D + 

A(D + 

A(D + 

A(D + 

A(Df 

A(Df 

A(Df 

A(Df 


vr+vr-) = -V cd V ud (T + E q ) 

0 On _ VcdVud . 


7T 7T = 


V2 


C + E q ) 


A) = + Cf q ) cos </> - V CS V U a cf a ^ - V cd V ud E q cos <j> 

A) = ^r^(- A + C# ) sin 0 + F CS K A ^ Sin 0 

w) = Vcd J^ d ( c ( q + A cos2 0 + VcsVus(-C( a + Ap, sin 2 </>) 


W) = VcdVud{E t 


sin2 ^> , (r f + r f t sin A + y y r_r/ + cJ _ 

<? 2 4 J + A ^2 + A ^2 


cos - 


— E s sin 2(f>) 


K+K~) = V CS V US (T + E q ) 

K°K°) = -{V cs V us E q + V cd V ud E s ) 

> 7 r+ 7 r o ) = _^Kd (T + c) 

v2 

’ A) = Vcd J^ d (t£> + Cf q + 2A q ) cos (f) - V cs V us Cl s sin <f> 

' 7r + ^? / ) = ~ ,d ^ d (T%° + + 2^4 ? ) sin <f> + VcsVusC^ cos cf) 

► A' + j?°) = + K,K. s r 

► tt+A 0 ) = -{V cd V ud T + F CS KA) 

► 7r°A'+) = - -^(KdKdC - V cs V us A q ) 

► A) = (V cd V ud Cf lq + V cs V us A q )^t - V cs V us (tE> + Cl + A a ) sin0 

► A') = (VcoVudCl + V cs V us A q f-^t + V CS V US (T$ + + As) cos <f> 

The A>° —> 7T7r SCS decays obey the following triangular Isospin relation: 
A(D° —y tt + tt~) + \[2A(D° —> 7r°7r°) = \/2A(.D + —> 7r + 7r°) 


for all the n D —>■ PP decays, the FSI corrected amplitudes or the militarized’ amplitudes, 
A with % = 1, ...7i are given by 2 

N 

-41 = E 5 I /2 A» < 12 ) 

k=l 

2 For part of the formalism used in this section, we closely follow Refs. [66, 94] 
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CF Decays 


(D° -> K~tt + ) = -V cs V ud {T + E q ) 

(D° -+ AV) = -V cs vJ^M 

(• D° -> K° V ) = V cs V ud ((-C% -E q )^t + E s sin </>) 

(D° ->• W) = VcsV ud ((-C^ -E q f-^t- E s cos0) 

(D+ _> R\ + ) = -V cs V ud (T + C ) 

(Df -> A°A+) = -V cs V ud (C + A a ) 

(Df n + T]) = VcS J^ d (T^° sin (f>-Ag cos (f>) 

(D+ -> tt+t/) = (-T^° cos </> - A, sin </>) 

The -D 0 —>• /i 7r CF decays obey the following triangular Isospin relation: 

A(D° -> A"7t + ) + V2A{D° -> AV) = t 4(D+ -> A°vr + ) 

DCS Decays 

(D° ->■ A+tT) = V cd V us (T + E q ) 

(D° ->■ AV) = VcdVus ^^^ 

(D° -> K°r]) = V cd V us {(C% + E q )^t- E s sin0) 

(D° ->■ A 0 ?/) = V cd V us ((C^ +Eq ) S ^ J t+ E s cos </>) 

(D+ ->■ A<V+) = F cd F us (C + A,) 

(D+ -> A+tt 0 ) = V cd V u ^^ 

(D+ A+tj) = -V cd V us {(T^ +A q f-^t- sin0) 

(D+ ->■ A+r/) = -VcdVusdTp + + A s cos 0) 

h V2 

(D+ -> A+A°) = V cd V us {T + C ) 

The D° —>• I\tt DCS decays obey the following quadrilateral Isospin relation: 
A{D° -> A+tT) + V2A{D° -> K°ir°) = A{D + -> A 0 tt + ) + V2A(D + -> A + ^°) 
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where, S is the strong interaction matrix and k — 1, ...n,n + 1, ....N, stands for all possible 
states that can rescatter into the PP states. In the heavy quark limit the hard rescattering 
dominates, in which case, the sum can be interpreted to be over all intermediate states of 
partons and the number of these states will hence be very large. Parton hadron duality will 
then permit this estimation. These corrections are incorporated into the hard scattering 
contributions in the QCD factorization approach of Ref. [73] for the case of B meson decays. 
For the case of charmed meson decays, since m c is not large enough, we include these NF 
corrections in the modified Wilson coefficients, in terms of parameters Xi and Xi- However, 
some residual long distance FSI’s may be left which are particularly important for charmed 
meson decays, due to the nearby resonances. This residual rescattering is considered in the 
limited set of D —> PP decays, to which the duality cannot be applied and therefore these 
effects may not be incorporated in the NF corrections. 

The S matrix in Eq. (12) can be written in terms of a residual matrix (<S res ) for the 
rescattering among the D —> PP states alone and the scattering matrix which accounts for 
the hard rescattering from all possible hadronic states into these channels, resulting in the 
factorization amplitudes as, 


S ik = y>i)ij(<S 2 )jk, where, Si = S res and 
j=i 

N 

S 2 — SpiS.with A**' — W '-Sp i,: .-I, resulting in 

•4F = x>;/ s 2 Mp 


k=i 


k=i 


Any S matrix can be written in terms of a real, symmetric K matrix as, S = 
Hence, the unitarized amplitudes in Eq. (13), may be written as, 


(13) 

(1 - iK)“‘(l + iK). 


^ = £«1 - iK)" 1 ),^ . (14) 

k=l 

The K matrix parametrization has the advantage that the resonances coupling two body 
channels are represented by poles in the K matrix. The summation in Eq. (14), corresponds 
to summing the geometric series, where the final state hadrons are produced from scattering 
via resonance at different orders, starting from zero, i.e., directly from the decaying meson 
without the resonance contribution, resonant rescattering occurring once, twice and so on. 
While such coupled FSI’s has been considered in the past in many papers [31, 64, 71, 76, 95] 
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most of the papers on charm decays further assume SU(3) to relate the parameters of the 
coupling matrix. Since SU(3) is broken, we prefer to use the measured decay rates of the 
resonances to various channels to fix the K matrix parameters as far as possible and the ones 
not measured are left as parameters to be determined by fits of all the theoretical branching 
ratios to the observed values. 

For each of the SCS, CF as well as DCS modes, states with the same isospin are coupled 
together. In general the K matrix coupling three channels will have the form: 

ill'll \fk1k2Y12 yfk\k 3 r 13 
\fk2klY2i k 2^22 a/W 23 

VhhTn 32 A; 3 r 3 3 , 

where, nines denotes the mass of the resonance through which the different channels are cou¬ 
pled and hi, k 2 and /c 3 are the cm momenta of the 3 decay modes. There are six independent 
parameters . To reduce the independent parameters to a manageable number, we impose 
the requirement that the diagonal cofactors of K(s) vanish (or equivalently, detiC(s) = 0 ). 
This leads us to three conditions, 

t 2 12 = r u r 2 2, y 2 13 = r n r 33 and r2 3 = r22r 33 . (15) 

The T' -s are related to the the partial decay width of the resonance to the i th channel. 

To illustrate this, we consider first the case of isospin zero states of SCS decay modes of 
D° meson. The isospin zero combination of the 7r + 7T _ and 7r°7r°, K + K~ and K°K°, and 
the 7777 modes, are coupled via Eq. (14) with the /o(1710) pole in the K matrix. Hence, for 
this specific case of coupling of the isospin zero states, in the K matrix, nines — 1-720 GeV, 
k\ = ^(rn D o 2 — 4m?o ), k 2 = |( m D o 2 — 4ni 2 J<(i ), k 3 = k(m D o 2 — 4m^) and we have 

r(/ 0 ^vrvr) = ^, r(/ 0 ->■ KK) = and r(/ 0 ^ 7777 ) = 

WlRes WlRes WlRes 

Experimentally only the two ratios of the decay rates, r(/ 0 —> KK)/F(f 0 —> rnr) and 
r(/ 0 —> KK)/T{f 0 —> rjr 7 ) have been determined. Hence we keep g pe = T(/ 0 —> KK) as a 
parameter, to be determined from fits of our theoretical estimates to the observed branching 
ratios. 

Similarly, for the 1=1 case, we take the Oo(1450) resonance with nin es = 1.474 GeV and 
r Res — 0.265 GeV, to be responsible for the rescattering among the channels KK, nr) and 


K(s) = 


\ m Res - S > 
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nrj', to which this resonance decays. Here again, the decay rate T(ao —» nrj) is not yet 
accurately measured and is treated as a parameter (h pe ) that may be predicted from the 
fits of all the branching ratios of the D —> PP modes to experimental data. Note that the 
KK , nr} and nr/ states appear as final states not only of SCS D° and D + decays, but also 
in the CF decays of the Df decays. For all these three sets of decays, the same K matrix 
(apart from tiny modifications in the cm momenta and the mass-squared of the decaying 
meson) will suffice, and more importantly with the same one unknown parameter, while 
many additional observables (all the branching ratios of these D°, D + and Df) will get 
added to the y 2 fit. If this one same parameter, along with the other unknowns in our 
analysis can simultaneously explain all the observed data, it would indicate that our naive 
technique of incorporating the FSI effects is satisfactory. 

We also couple the isospin 1/2 states of the Kn, Kr} and Kr/ channels, that are the final 
states in the SCS decays of D+, CF decays of D° and DCS decays of D° and D + mesons. 
Here we use the A'q( 1950) resonance with rrif> es = 1.945 GeV. Only the branching ratio, 
T(K* —>• Kn)/T tota i has been measured. We take the other two decay rates, T(K/ —$■ Kr /) 
and r(/lg —> Kr f) as parameters (j pei and j pe . 2 ) that can be determined by the overall fits 
of all the branching ratios to data. 


IV. NUMERICAL ANALYSIS AND RESULTS 

To estimate all the possible sets of coupled channels, the isospin decomposition of all 
the SCS, CF and DCS modes are listed on the next page. Here A ( ' U} denote the bare or 
un-unitarized (unitarized or FSI corrected) amplitudes respectively, for each of the decay 
modes, while A/ lode(u > denotes the corresponding ununitarized (unitarized) isospin, I = i 
amplitudes for those modes. With absence of resonances in particular isopin components 
with the right quantum numbers in the vicinity of the charmed meson masses, some of the 
isospin components of many modes remain un-unitarized. 


With all the unitarized isospin amplitudes, we construct the corresponding unitarized 
decay amplitudes for all the decay modes. The decay rates for all the D —> P 1 P 2 are then 
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SCS Decays 


AW(D° - 

> 7T+7T-) = VZA$* + V2A™ {U) 

A {U) (D° - 

> 7r°7r°) = 2 AY - Al n(U) 

A {U) (D° - 

> vr °n) = V3Ar (U) 

a( u \d° - 

■> 7r°r/) = y/3Ai'' l ' b ' 

A {U) (D° - 

■> ??r/) = V3Aq' i< ' U ' > 

A {U) (D° - 

> W) = Va4^ 

A {U) (D° - 

>K+K-)=yJl(A? Km + A™ m 

A {U) (D° - 

t K°K°) = J\ (A« Km - 4 K(U> ) 

A i - U \D + - 

-> ^ + 7T°) = 3^ 

A { - u \d + - 

■+ A^ + A 0 ) ^ Af^ 

A<M(D+ - 

-+ 7 r+7 ? ) = ,?(C/) 

a { - u \d + - 

-* TT+t/) = ^ +??,(C/) 

A (C/) (A>+ - 

^ +K °^TAf + ^l A T u) 

A (C/) (A>+ - 

- n°K + ) . y|^f - <"> 

A (C/) (A>+ - 

■* K + rj) S ^ + ’' (t ' ) 

2 

A (C/) (A>+ - 

■* *V) = 


2 


CF Decays 


AV\d° - 

■> A" 7T + 

) = \a^ + \a^ { u) 
0 2 0 2 

AW(D° - 

+ A"°7T 0 ) 

0 2 2 

A (C/ )(D° - 

+ A^r?) = 

- 

A iu) (D° - 

+ K°rj') ■ 

- 
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A( u \d + -> K°n + ) = A§ n+ 

2 

H (C/) (D+ -> A°A^ + ) = A^ K{JJ) 
AM(D+ ^ir + r l ) = Af 1,{U) 

aW( D + ^ ttV) =-4T V(C7) 



DCS Decays 


aM(d° - 

» A+7T-) = ^ A\ 

O 2 

V2 A Kn(U) 

Vs '§ 

A<W(D° - 

» K°ti r°) = + ■ 

O 2 

1 A Kn(U) 

Vs h 

AV\D° - 

> K°n) = a k V { u) 

2 


A<W(D° - 

> K°r/) = Af' m) 

2 


AV\d + - 

4 A-V+) = 

3 2 \/3 2 

a ( - u \d + - 

-* K + 7T°) = \A^ + - 
6 2 

1 .K*+(U) 

Vs h 

AV\d + - 

k :+,) = 

2 


a ( - u \d + - 

4 K+rf) = 

2 


AW(D+ - 

+ k + k°) = A[a? k + a$ k ) 


calculated as, 

r(D PxP 2 ) = PiP 2 )| 2 . ( 16 ) 

Here p c is the centre of mass momentum of the mesons in the final state given by 

The theoretical branching ratios for each of the decay modes of the D°, D + or the Df 
mesons are then obtained by dividing the corresponding decay rates by the total decay 
widths of these mesons. We then perform a y 2 fit of these theoretical branching ratios with 
the experimentally measured branching fractions, estimating all the unknown parameters 
from the best fit to data. 
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The unknown parameters in onr study are: the four parameters representing the NF 
corrections, yi, X 2 and their respective phases (f) 1 , d 2 , four parameters: Xq, s an d Xq )S depicting 
the strength of the W-exchange and W-annihilation amplitudes with distinct strengths for qq 
and ss pair production, one unknown in each of the isospin zero and isospin one K matrices 
coupling modes from decays of D° and D + mesons, two parameters in the isospin half K 
matrix coupling various decay modes of Df, one parameter A, representing the momentum 
of the soft degrees of freedom in the charmed mesons, that is used to define the scale for 
each of the individual decay modes, making a total of 13 unknown parameters. On the 
other hand out of all the 33 decay modes considered, 28 have been measured, resulting in 
sufficient observables to deterimine all the unknown parameters and give predictions for five 
of the branching fractions of DCS modes that are not yet measured. 

Apart from the experimental errors in the observed branching ratios, the calculated errors 
in our theoretical branching ratio estimates arise from the errors in the form factors, the 
77 and 7 / decay constants, the 77 — rj' mixing angle and the errors in the measured decay 
widths of the various resonances into the different channels that are included in our y 2 fits. 
Errors due to the other theoretical inputs, like meson masses, decay constants of pion, Kaon, 
charmed mesons, CKM mixing elements (involving the first two generations) are negligibly 
small. 

As discussed in Sec. IIC, the z-series expansion has been obtained for the form factors, 
keeping the first two terms of this series. The coefficients Oo and a\ of these two terms 
are functions of normalization and shape parameters /( 0 ) and /3, obtained from lattice 
results. The errors in these lattice parameters are used to obtain the errors in the expansion 
coefficient functions and then propagated to get the errors in the form factors. We find that 
the errors in the form factors vary from ^5% to ~25%. 

In the tables IV, V and VI, we list the values of the Branching Ratios of the all the 
SCS, CF and DCS D —> PP modes obtained from our analysis, after incorporating the FSI 
effects (shown in the 2 nd column), as well as in the absence of the FSI (column 3), absence 
of annihilation (column 4) along with the corresponding observed experimental branching 
ratios (column 5) given in PDG [91]. We also predict the B.R.’s for a few DCS modes that 
have not been experimentally measured yet and are given in the second column of Table VI 

After incorporating all the errors, the y 2 minimization results in the following best fit 
values of all the parameters: 
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TABLE III. Parameter best fit values 


Name 

Values 

Name 

Values 

Name 

Values 

A 

0.625645 

Jpe 1 

0.0000239368 

Y A 

A,q 

132.685 

Xi 

-2.68215 

Jpe 2 

0.096456 

y A 

As 

193.447 

X2 

2.23605 

Y E 

A,q 

-334.805 

<f>i 

0.302258 

9pe 

0.0471262 

y E 

As 

-81.3363 

<t>2 

2.87681 

hpe 

0.118834 






V. CONCLUSIONS 

For several decades, various ratios of decay rates of many of the D —y PP modes remained 
to be a puzzle as these were expected to be one in the 577(3) limit, but the measured values 
exhibited large deviations from unity. 

We have evaluated the bare amplitudes of all the D —> PP modes using factorization, 
however, we add non-factorizable corrections, weak annihilation and exchange contributions 
as parameters, and in the hadron matrix elements, the q 2 dependence of the form factors 
involved are evaluated using the z-expansion method and finally, resonant final state interac¬ 
tion effects are incorporated. The parameters of the K matrix coupling the various channels 
are defined using the measured decay widths of the resonances (where available) and those 
unobserved, are left as parameters to be fitted from all the measured 28 D —>■ PP branching 
ratios. Our best fit has a y 2 /degree of freedom of 2.25, which is an improvement over the 
previous results in Refs. [59, 75]. 

• We are able to get reasonable fits to almost all the observed branching ratios. In 
particular, our branching fractions for D —> KK , D — )■ tttt modes that have been a 
long standing puzzle are in agreement with the corresponding measured values. 

• We have evaluated the cosine of the strong phase difference between the unitarized 
amplitudes for D° —> K~n + and D° —> K + n~ and obtain, cos 8kk = 0.94 ± 0.027. 
This result is consistent with the recently measured BESIII result, cos 5f[n = 1.02 ± 
0.11 ±0.06 ±0.01 [96], 

• The mode D° —> K°K° does not have any tree or colour suppressed contributions, but 
can come only from W-exchange. In fact, there are two exchange contributions, one 
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TABLE IV. D —> PP SCS B.R.’s, Columns 2 and 3 show our results with annihilation included, 
for the with and without FSI cases respectively, while column 4 which displays results without 
annihilation, includes FSI. 


Modes 

With FSI 

Without FSI 

Without Ann 

Experimental Value 

D° -> TT + TT~ 

(1.44 ±0.027) x 10“ 3 

(4.35 ± 1.67) x 10 -3 

(4.02±1.75x) x 10 -3 

(1.402±0.026) x 10 -3 

D° -> 7r°7r° 

(1.14 ±0.56) x 10" 3 

(3.66 ± 1.43) x 10 -3 

(2.04 ±0.79) x 10" 3 

(8.209 ±0.35) x 10“ 4 

D° -> K+K~ 

(4.06 ±0.77) x 10“ 3 

(4.27 ±2.34) x 10“ 3 

(6.78 ±3.08) x 10“ 3 

(3.96 ±0.08) x 10“ 3 

D° -> K°K° 

(3.42 ±0.52) x 10 -4 

(5.61 ±0.00) x 10 -4 

(2.80 ±0.84) x 10 -4 

(3.4 ±0.8) x 10" 4 

D° —> ir°r) 

(1.47 ±0.90) x 10 -3 

(6.47 ±2.98 x 10~ 3 

(3.25 ± 1.51 x 10 -3 

(6.8 ±0.7) x 10- 4 

D° —> Tr°r)' 

(2.17 ±0.65) x 10 -3 

(3.81 ± 1.43) x 10 -3 

(1.85 ±0.79) x 10 -3 

(9.0 ±1.4) x 10 -4 

D° —> 

(1.27 ±0.27) x 10 -3 

(1.32 ±0.41) x 10 -3 

(1.34 ±0.29) x 10 -3 

(1.67 ±0.20) x 10 -3 

D° —> i]r)' 

(9.53 ± 1.83) x 10 -4 

(1.04 ±0.27) x 10 -3 

(5.38 ± 1.63) x 10 -4 

(1.05 ±0.26) x 10 -3 

D + 7T + 7T 0 

(8.89 ±4.51) x 10 -4 

(8.70 ±6.70) x 10 -4 

(9.73 ±3.94) x 10 -4 

(1.19 ±0.06) x 10 -3 

D + -> K+K° 

(3.75 ± 0.63) x 10 -3 

(1.02 ±0.37) x 10" 2 

(1.99 ±0.56) x 10 -2 

(5.66 ±0.32) x 10 -3 

D + 7 T + r] 

(4.72 ±0.21) x 10 -3 

(2.34 ± 1.26) x 10 -2 

(1.66 ±0.77) x 10 -2 

(3.53 ±0.21) x 10 -3 

D + —>• 7 T + rj' 

(6.76 ±2.19) x 10 -3 

(3.00 ±0.76) x 10" 2 

(9.78 ±3.35) x 10 -3 

(4.67 ±0.29) x 10 -3 

D+ -f vr+A' 0 

(1.96 ±0.90) x 10" 3 

(1.46 ± 1.10) x 10 -3 

(1.32 ± 1.01) x 10 -3 

(2.42 ±0.12) x 10 -3 

D+ -> 7t°A" + 

(8.17 ±4.64) x 10“ 4 

(1.74 ± 1.00) x 10" 4 

(1.01 ±0.54) x 10" 3 

(6.3 ±2.1) x 10“ 4 

D+ K+r) 

(1.50 ±0.75) x 10" 3 

(6.40 ±4.52) x 10" 3 

(2.23 ± 1.82) x 10" 3 

(1.76 ±0.35) x 10" 3 

D+ K+r] 1 

(7.07 ±0.49) x 10“ 4 

(2.09 ± 0.87) x 10 -3 

(0.57 ±0.47) x 10" 4 

(1.8 ±0.6) x 10 -3 


appearing with a dd and the other with an ss, which under exact SU(3) symmetry 
would cancel each other, resulting in a null amplitude. However, since our parameters 
for these two contributions are distinct, our bare amplitude for this mode is small but 
non-vanishing. There have been speculations [72, 76] that this mode can arise just 
from final state interactions, even in the absence of a weak exchange contribution. 
However, from Table IV it is clear that without the exchange contribution we are 
unable to generate a large enough rate: both final state interaction and the exchange 
contribution are necessary for consistency with the measured branching fraction. 

• We have also evaluated the four ratios of amplitudes that had been specified in a 
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TABLE V. D —> PP CF B.R.’s, inclusion of Annihilation/FSI in our branching ratio estimates 
shown in various columns is the same as specified for Table IV. 


Modes 

With FSI 

Without FSI 

Without Ann 

Experimental Value 

D° -> K~tt+ 

(3.70 ± 1.33) x 10“ 2 

(8.83 ±2.47) x 10“ 2 

(5.63 ± 1.81) x 10“ 2 

(3.88 ±0.05) x 10“ 2 

D° -» K°ir° 

(1.88 ±0.99) x 10“ 2 

(1.29 ±0.44) x 10 _1 

(3.30 ± 1.47) x 10 -2 

(2.38 ±0.08) x 10“ 2 

D° -> K°r] 

(1.59 ±0.48) x 10“ 2 

(0.97 ±0.33) x 10" 2 

(1.09 ±0.34) x 10“ 2 

(0.958 ±0.06) x 10“ 2 

D° -> K°t]' 

(2.29 ±0.43) x 10“ 2 

(2.06 ±0.30) x 10 -2 

(2.45 ± 0.47) x 10“ 2 

(1.88 ±0.1) x 10“ 2 

D + -> K°tt+ 

(3.42 ± 1.78) x 10“ 2 

(1.35 ± 1.12) x 10 _1 

(5.25 ± 3.34) x 10“ 2 

(2.94 ±0.14) x 10" 2 

D+ -> K°K+ 

(5.65 ± 1.29) x 10 -2 

(1.70 ±0.79) x lO” 1 

(1.35 ±0.53) x 10” 1 

(2.95 ±0.14) x 10“ 2 

Df — > 7 T + r] 

(2.26 ±0.82) x 10“ 2 

(0.78 ±0.56) x 10“ 2 

(2.14 ±0.90) x 10“ 2 

(1.69 ±0.10) x 10 -2 

Df —» 7 T + r]' 

(2.64 ±0.78) x 10“ 2 

(3.73 ± 1.52) x 10“ 2 

(2.52 ±0.85) x 10“ 2 

(3.94 ±0.25) x 10“ 2 


TABLE VI. D —> PP DCS B.R.’s, inclusion of Annihilation/FSI in our branching ratio estimates 
shown in various columns is the same as specified for Table IV. 


Modes 

With FSI 

Without FSI 

Without Ann 

Experimental Value 

D° -> K + tt~ 

(1.77 ±0.88) x 10“ 4 

(3.71 ± 1.33) x 10 -4 

(2.48 ± 1.07) x 10“ 4 

(1.38 ±0.028) x 10“ 4 

D° -> A" 0 7t° 

(2.11 ±0.26) x 10“ 4 

(3.70 ± 1.35) x 10“ 4 

(0.68 ±0.46) x 10“ 4 

- 

D° -> K°r] 

(0.94 ±0.45) x 10“ 4 

(0.28 ±0.10) x 10“ 4 

(0.96 ±0.32) x 10“ 4 

- 

D° -> K 0 ri' 

(8.02 ±3.32) x 10“ 4 

(0.59 ±0.08) x 10“ 4 

(9.22 ± 1.61) x 10“ 4 

- 

D + -f K°n + 

(3.27 ± 1.86) x 10 -4 

(1.19 ±0.55) x 10" 3 

(3.51 ±2.11) x 10“ 4 

- 

D + -f K + n° 

(3.07 ± 1.02) x 10“ 4 

(2.15 ± 1.17) x 10“ 4 

(3.27 ± 1.39) x 10 -4 

(1.83 ±0.26) x 10 -4 

D + K + r] 

(0.98 ±0.26) x 10“ 4 

(1.04 ±0.23) x 10“ 4 

(0.89 ±0.27) x 10“ 4 

(1.08 ±0.17) x 10 -4 

D+ -> K+ri' 

(1.40 ±0.39) x 10 -4 

(1.82 ±0.18) x 10“ 4 

(1.35 ±0.39) x 10“ 4 

(1.76 ±0.22 x 10' 4 ) 

D+ -> K+K° 

(7.84 ±2.31) x 10“ 4 

(0.68 ±0.09) x 10“ 4 

(0.72 ± 0.44) x 10” 4 

- 


recent paper [97]. In SU(3) limit these are all expected to be unity. Our theoretical 
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estimates for these ratios are given below: 


Ri = 
Ra = 
Rs = 


\A(D° -2 K + n-)\ 

_ 1 _1_ Li _— i 07 4- n qo 

\A(D° —>• n + K~ | tan 2 6 C ’ 

I A(D° -2 K+K~)\ 

i_1_hi — 1 27 ± 0 42 

|7L(A° ^7T+7T-| ’ 

|7L(A>° -> Jl + A'-)| + |7L(A>° -2 7T+7r-| 

171(A) 0 -9 7T+/1 i tan0 c + |7L(A)° -9 K +7T")| tan" 1 0, 


A 4 = 


/171(A) 0 -2 AT+AT-) | \A(D°n+n- 


\A(D° -2 7r+A'-||7l(A 0 -2 AT+tt-)| 


= 1.19 ±0.26. 


1.19 ±0.28, 


Furthermore, the following combination of these ratios is expected to be vanishing up 
to 4 th order in U-spin breaking, 


AA = R -3 — Ra ± — 


V 


V^Ri 


Using our unitarized amplitudes, we hnd the central value of A R to be indeed very 
tiny, however, with a large error. 


A R= -0.000013 ±0.006. 


Our theoretical errors (in form factors, K-matrix parameters, etc.) are propagated to 
evaluate the errors in the ratios Ah, Ah, A 3 , A 4 and finally in A A. 

• We would like to mention that in many modes involving and 7/, we have additional 
terms in our amplitude due to our distinction of the different form factors, compared 
to, for eg., those that appear in Ref. [59]. A naive look at the colour suppressed 
diagrams for D —> ir°r}(r)') will indicate that the contributions from the case where 
the spectator is part of the 7r°, and that, where it constitutes the 77 ( 7 /) must cancel. 
However, in terms of the specific decay constants and form factors, one is proportional 
to —f n F® nq (m^) while the other is proportional to f Vq F^ n (m^), which are unequal 
and hence must survive. 

• While the Particle Data Group [91] does not include a world average for r(/ 0 (1710) —> 

KK) but it does list two values for the ratio : 0.36 ±0.12 (Ref. [98]) and 

0.381^9 (Ref. [99]); our fit value of g pe (= r(/o(1710) — > KK)) corresponds to 0.35 
for the branching ratio, which is consistent with these values. 
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• Our theoretical errors are rather large and could be reduced in future with more precise 
form factors available either from measurement of scmileptonic D, D s modes at BES 
III, where if even the lepton mass could be incorporated (eg. by looking at modes with 
muon in the final state), then the q 2 dependent F 0 could be known, or with improved 
lattice studies, specially for D —>• rj', D s —> K and D s —> 77,7/. 

• Accurate measurements of the decay widths of the resonances (used for the final state 
interactions) to many of the coupled channels can reduce the theoretical uncertainties 
and possibly allow for better fits to data. For example our fits seem to indicate a 
rather large value for h pe or the width of ao(1450) —y n°r). This seems to result in 
larger branching fractions for many of the isospin one modes. Future measurement of 
this width can help reduce this uncertainty and perhaps result in better fits to data 
for these modes. 

• Out of the 28 observed PP modes, we are unable to fit 7 of the modes well. Many 
of these modes involve rj or rf in the final state. Including a gluonium component 
in the 77, rf states may possibly be one way of improving these fits. This, along 
with improved form factor measurements, observation of decay rates of the resonances 
(playing a role in final state interactions) to these decay modes, as mentioned in the 
last two points above, could go a long way in improving our fits. One glaring misfit 
is the mode D + —y K + K°. This mode does couple to 7r + ?7,7r + // modes and hence, 
may possibly improve, along with the improvements in those modes. Finally, perhaps 
a more sophisticated statistical analysis may also play a role in improving our results, 
which we hope to carry out in future [100]. 


Appendix A: The Wilson Coefficients 

The Wilson coefficients used in the evaluation of all the bare amplitudes have been 
calculated at the final state hadronic scale. This allows for an additional SU( 3) breaking 
effect and they have been incorporated using the procedure outlined in Ref. [75]. The Wilson 
coefficients at lower scale are calculated in the Ref. [101]. The essential steps are following. 

• In the first step, the Wilson coefficients Ci(mw) at weak scale are calculated by re¬ 
quiring the equality of the effective theory with five active flavors q — u,d, s, c, b onto 
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the full theory. 


• Next, the coefficients undergo the evolution from the scale rriw to // through the 
equation 


C(fj) = U 5 (fi,m w )C(m w ). (Al) 

• In the next step, coefficients are calculated at the scale of b quark 


C(m b ) ->• Z(m b ). (A2) 

• Now, the Wilson coefficients can be evaluated at required scale (/i c or fihadron ) through 
the equation 


C(fj) = U 4 (n,m b )Z(m b ). (A3) 

In the above steps, U 5 and f/ 4 are the 2x2 and 7x7 evolution matrices for five 
and four active flavors respectively. The Z(m b ) is given in the Eqs.A.7 to A.10 of the 
Ref. [101]. 

The explicit expressions of the Wilson coefficients obtained after following the above steps 
are given in the Ref. [75]. They are: 


C\(n) = - 0.2334a L444 + 0.0459a a7778 + 1.313a a4444 - 0.3041a-°' 2222 , 
C 2 (n) = 0.2334a 1 ' 444 + 0.0459a a7778 - 1.3 1 3a a4444 + 0.3041a-°- 2222 


in terms of the running coupling constant a: 

/3ilnln([i 2 / A^jg) 

A|^) 

with the coefficients 


a = Oi s (n) = 


An 


(3 0 ln(n 2 / A| 


MS' 


fa 


33-2/ 

3 


ft = 102 - 


We take active flavour number Nf = 3, and the QCD scale A jjg = = 375 MeV. 

Again, note that the scale dependent strong coupling constant a s (/x) is evaluated at the final 
state hadronic scales for each individual decay, to take care of the SU(3) breaking. 
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Appendix B: Series Expansion Method for Form Factors 


Using the analytic properties of F(q 2 ), a transformation of variable is made which maps 
the cut on the q 2 plane onto a unit circle \z\ < 1, where 


z{tt 0 ) 


y/t+ — t — y/t+ — to 
\/t+ t + \Jt + — to 



where, 


t± = (mo ± mpj 2 and t 0 = t + 1 



This transformation allows the form factors to be given by an expansion about q 2 = to, given 
as 

^ inf 

m = pmt,t o) Z> (f oMM')‘, 

given also as Eq.( 8) in the text. The function P{t) in the above is 1 for D —> n form factors. 
For D s —y r/ and D —> K form factors, P(t) = z(t,Mj ) (where Mp* ) is the nearest 0 + 
resonance mass). The outer function cj) is given by [102] 


0(t,t o ) = 


3t+t- fz(q 2 ,0)\ 2 fz(q 2 ,t 0 )\ 1/2 fz(q 2 : t-)\ 1/4 


32vry 0 


to - q 2 


t_ — g 2 


where yo has been calculated [102] using OPE and is given by: 

1 + 0.751 a s (m c ) 


Xo = 


8tt 2 


y«+ - <i 2 


For simplicity, we ignore condensate contribution which is of the order 0{m c 3 ) and 0(;m c 4 ) 
The strong coupling at charm scale is computed with the package RunDec [103]. 
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